ellipsenm: ecological niche model
This is the final report produced during the execution of the function ellipsoid_model from the ellipsenm package. Following, a complete description of the ecological niche modeling process for the species pmorio is presented. An ecological niche, from a Grinnellian perspective, is the set of environmental conditions that allow a species to maintain populations for long periods of time, without immigration events (Peterson et al. 2011). Models created with the ellipsenm package are ellipsoid envelope models and assume that a species ecological niche is convex, has an only optimum, and the species response to each variable covaries with the response to other variables (see Brown 1984, Jimenéz et al. 2019, Osorio-Olvera et al. 2019). Mahalanobis distances are used to represent how far is each combination of environmental conditions from the optimum (ellipsoid centroid). Suitability values result by default from a multivariate normal transformation of Mahalanobis distances. Therefore, maximum values of suitability will be close to the centroid and minimum values will be close to the border of the ellipsoid. To obtain binary layers from the ones written in the output_directory for suitability predictions, reclassify the values in the layers as follows: 0 = 0 (non-suitable, outside ellipsoid), and values other than 0 = 1 (suitable, inside ellipsoid).
Note 1: If only one replicate was performed, results are obtained using the hole set of data. If multiple replicates were performed, results are produced per replicate; as each replicate is build with a sub-sample of the data, results may differ among replicates. For replicated models, mean results are shown in this report, results for all replicates, as well as for mean, minimum, and maximum ellipsoidal models can be found in the output directory defined in the function ellipsoid_model. If the model was performed with replicates, further raster calculations can be performed to show model variability (e.g., suitability range or other indicators of model variability; see Peterson et al. 2018)
Note 2: This is the report for a model projected into scenarios other than the one used for calibration. Therefore, other results correspondent to predictions for each of the projection scenarios were created in the output directory defined in the function ellipsoid_model. This report includes graphical and numeric results for only one of these projections (generally the first of all projections if multiple, or the only projection if a single one was performed). The reason for presenting only one projection is the prevention of RAM-related issues.
Note 3: If more than tree dimensions were used to create the model, up to three of them will be used in graphical representations.
Tip: To see two dimensional figures in their actual size click on them. Three dimensional figures are interactive.
Summary
The model presented in this report was produced for the species pmorio to characterize its ecological niche using an approach based in ellipsoid envelopes. The modeling process was performed with 1 replicate. Based on the level used to produce the ellipsoids (99) the error (E) assumed for the occurrence data was 1%. This is, 1% of the occurrence data is considered as potentially erroneous (i.e., environmental outliers); therefore, this data is not included in the ellipsoid envelope model for the species ecological niche.
This model was projected to 1 scenario(s). The results for all projections can be found in the output directory defined in the function ellipsoid_model under the name(s): projection. Graphical and numerical results are presented only for the projection scenario.
Data used for modeling
A total of 108 occurrences and 6 variables (wc2.1_2.5m_bio_2, wc2.1_2.5m_bio_3, wc2.1_2.5m_bio_5, wc2.1_2.5m_bio_13, wc2.1_2.5m_bio_15, and wc2.1_2.5m_bio_18) were used to produce the models.
Data in geographic space
The plot below shows the geographic arrangement of occurrence data on the variable wc2.1_2.5m_bio_2. Colors correspond with values of the variable according to figure legend.
Figure 1: Representation of species occurrence data in geographic space. Occurrence data is represented in black. Grey lines represent country borders.
Data in environmental space
Calibration area
The plot below shows the arrangement of the occurrence data in environmental space considering 3 dimensions for the area of calibration.
Figure 2: Representation of species occurrence data in environmental space for the calibration area. Grey points respresent environmental conditios in the calibration area; blue points respresent occurrence data.
Representation of species occurrence data in environmental space for the calibration area. Grey points respresent environmental conditios in the calibration area; blue points respresent occurrence data.
Projection scenario
Below a plot is showing the arrangement of the occurrence data in environmental space considering 3 dimensions for the projection scenario.
Figure 3: Representation of species occurrence data in environmental space for a scenario of projection. Grey points respresent environmental conditios in the calibration area; blue points respresent occurrence data.
Representation of species occurrence data in environmental space for a scenario of projection. Grey points respresent environmental conditios in the calibration area; blue points respresent occurrence data.
Numeric results
Numeric results of the ellipsoid envelope model for pmorio are presented below.
Ellipsoid characteristics
The table below describes the characteristics of the ellipsoidal ecological niche model. The complete report of ellipsoids’ characteristics (e.g., centroid, covariance matrix, semi-axes length, etc.) can be found in the txt files referred in the table.
| ellipsoid_model | |
|---|---|
| Method | covmat |
| Level | 99 |
| Volume | 43915229328.51 |
| Other_metadata | ellipsoid_metadata.txt |
Prevalence in calibration area
Prevalence in environmental and geographic space are summarized below. Prevalence is defined as the proportion of space identified to be suitable for the species. The difference between the two types of prevalence is that in geographic space proportions are calculated using all pixes from the raster; whereas, in environmental space, only pixels that have distinct combinations of all variable’s values are considered (see Hutchinson’s duality in Colwell and Rangel 2009).
| ellipsoid_model | |
|---|---|
| prevalence_E_space | 0.7575385 |
| prevalence_G_space | 0.7575385 |
Prevalence in projection scenario
Prevalence for the projection scenario are presented below (if more than one projection was performed, prevalence in other scenarios can be found in the output directory).
| ellipsoid_model | |
|---|---|
| prevalence_E_space | 0.3572813 |
| prevalence_G_space | 0.3572813 |
Predictions in environmental space
Given the data and the characteristics that define the species ellipsoidal niche (see Data used for modeling and Ellipsoid characteristics) the results in environmental space are presented bellow.
Calibration area
The calibration area is represented by the variables given to fit the ellipsoid using the ellipsoid_model function. This area is recommended to be a region that has been accessible to the species for a relevant period of time (see Barve et al. 2011). Although most methods used in the ellipsenm package are not sensible to a background (environmental conditions in the calibration area), having a well defined accessible area helps with interpretations. Bellow, a 3 dimensional representation of the ellipsoidal niche model and the suitability values is shown in environmental space for the calibration area.
Figure 4: Representation of suitability values in environmental space for the calibration area. Occurrence data is represented in blue. Cold colors represent low suitabilities and warm colors high suitabilities.
Representation of suitability values in environmental space for the calibration area. Occurrence data is represented in blue. Cold colors represent low suitabilities and warm colors high suitabilities.
Projection scenario
A 3 dimensional representation of the ellipsoidal niche model and the suitability values is shown in environmental space for the projection scenario.
Figure 5: Representation of suitability values in environmental space for a scenario of projection. Occurrence data is represented in blue. Cold colors represent low suitabilities and warm colors high suitabilities.
Representation of suitability values in environmental space for a scenario of projection. Occurrence data is represented in blue. Cold colors represent low suitabilities and warm colors high suitabilities.
Predictions in geographic space
Given the data and the characteristics that define the species ellipsoidal niche (see Data used for modeling and Ellipsoid characteristics) the results in geographic space are presented bellow. Values of suitability can range from 0 to 1; when the maximum is not 1, it means that there is not environmental conditions equal to the ones at the centroid of the ellipsoid. Another reason why the range of suitability can change is that a different function of suitability was used in the ellipsoid_model function. Values of zero (coldest areas) represent places with environmental values outside the ellipsoid, values other than zero represent levels of suitability.
Calibration area
The geographic prediction in the calibration area is presented below. Colors correspond with values of suitability according to figure legend.
Figure 6: Representation of suitability values in geographic space. Cold colors represent low suitabilities and warm colors high suitabilities. Grey lines represent country borders.
Projection scenario
A geographic prediction in a projection scenario is shown below. This example corresponds to the scenario projection.
Figure 7: Representation of suitability values in geography for a projection scenario. Cold colors represent low suitabilities and warm colors high suitabilities. Grey lines represent country borders.
This report was produced using the ellipsenm R package; to cite R and this package use the following code:
The code for producing this report (if needed) can be found in the file enm_suitability_pr_report.Rmd; the data used is in enm_report_data.RData.
References
- Barve, N., Barve, V., Jiménez-Valverde, A., Lira-Noriega, A., Maher, S.P., Peterson, A.T., Soberón, J., Villalobos, F., 2011. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol. Modell. 222, 1810–1819.
- Brown, J.H., 1984. On the relationship between abundance and distribution of a species. Am. Nat. 124, 255–279.
- Colwell, R.K., Rangel, T.F., 2009. Hutchinson’s duality: The once and future niche. PNAS 106, 19651–19658.
- Jiménez, L., Soberón, J., Christen, J.A., Soto, D., 2019. On the problem of modeling a fundamental niche from occurrence data. Ecol. Modell. 397, 74–83.
- Osorio‐Olvera, L., Soberón, J., Falconi, M., 2019. On population abundance and niche structure. Ecography 42. https://doi.org/10.1111/ecog.04442.
- Peterson, A.T., Cobos, M.E., Jiménez‐García, D., 2018. Major challenges for correlational ecological niche model projections to future climate conditions. Ann. N. Y. Acad. Sci. 1429, 66–77.
- Peterson, A.T., Soberón, J., Pearson, R.G., Anderson, R.P., Martínez-Meyer, E., Nakamura, M., Araújo, M.B., 2011. Ecological Niches and Geographic Distributions. Princeton University Press, Princeton.